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O Abstract. The nilpotent fonnalism for the Dirac equation, outlined in previous papers, 

^ is applied to QED. It is shown that what is usually described as 'renormalization' is 

p I effectively a statement of the fact that the nilpotent formulation is automatically 

^ second quantized and constrains the quantum field into producing finite values for 

fundamental physical quantities. 

^ 1 Introduction 

ON 

^ In earlier papers ^"^ we have developed a formalism for the Dirac equation, based on a 

O nilpotent algebra, which has been applied in such cases as the hydrogen atom and the 

^ strong interacting potential. The idea has been to generate structures which explain 

^ various aspects of particle physics in a more fundamental way. Certain aspects of this 

^ algebra have also suggested that it might lead to a more fundamental understanding of 

the existence of apparent infinities in the QED formalism and their removal by the 
^ method of renormalization. It seems possible that a formalism which doesn't appear 

to require normalization, might not require renormalization either.^ Also certain 
> aspects of the algebra indicate that it has a built-in second quantization and a natural 

^ supersymmetry without requiring extra supersymmetric particles, which suggests that 

^ infinite sums must add up to finite values and that any necessary cancellations might 

be automatic. Our aim, therefore, has been to develop our formalism to incorporate 
the existing techniques of QED, and to discover if there are any aspects of the 
formalism which help either to explain or to remove apparent anomalies in the 
existing structures. A brief discussion of significant aspects of the nilpotent algebra 
will be followed by an outline derivation of those aspects of QED which are relevant 
to our discussion of the basis of the apparent divergences in the theory and their 
resolution through renormalization. 

2 The nilpotent formulation of relativistic quantum mechanics 

The nilpotent formulation requires the use of the 32-part algebra generated by 
combining the 4-vector units {i, i, j, k) with the quaternion units (1, i,j, k), where the 
two sets of terms correspond to the units of the fundamental parameters space-time 
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and mass-charge. The vector units (i, j, k), however, have the multivariate or 
quaternion-like properties of Pauli matrices. In effect, this means defining a 'full 
product' for two vectors a and b of the form 



ab = a.b + j a X b , 

and using the following rules for the multiplication of the vector units: 

i^ = f = k^=l 

ij = -ji = «k 
jk = -kj = a 
ki = -ik = rj . 

To construct the Dirac equation using this algebra, we begin with the classical 
relativistic energy conservation equation, 

—p^ — m^ = , 

which we factorize and supply with the exponential term e'^^' ' ^'^\ so that 

(± kE ± a p + ij m) (± kE ± ii p + ij m) e'^^' " p-"^ = . 

Then replacing E and p in the first bracket with the quantum operators, id I dt and 
-jF, to give / 3 A 

± ikr. + i V+ ijm \i±kE± ii p + ij m) e'^'^' ' p-"^ = , 

we obtain a quantum mechanical equation of the form 

± ± i V+ ijm y/ = , 



V 



where the wavefunction (or, in our terminology, the quaternion state vector) for a free 
fermion becomes 

y/= (± kE ± ii p + ij m) e'^'^' " "-""^ , 

and the following correspondences are seen between the conventional gamma 
matrices and the terms of the nilpotent algebra: 

f = ik, y = ii, f = ij, ^ = ik, ^ = ij. 



For convenience, we generally arrange the four terms incorporated in y/in a column 
(or ket) vector, while the four terms in the differential operator take up equivalent 
positions in a row (or bra) vector. 
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The quaternion state vectors in this formulation have remarkable properties. They 
are nilpotents or square roots of 0, and they have the properties of quantum field 
operators, creation and annihilation operators, and even supersymmetry operators. 
Because the differential operators are, in effect, identical to the quaternion state 
vectors, the theory already has the second quantization required of quantum field 
theory, and it has been rigorously shown that they are identical to quantum field 
integrals acting on vacuum. ' Another important aspect of the nilpotent state vector is 
that it automatically incorporates both nonlocality (because a superposition of two 
identical wavefunctions is instantaneously impossible with a nilpotent structure) and 
time-delayed local action (because nilpotents are square roots of a relativistic 
equation). 

Many important results have already been achieved with this formalism, but, for 
the purposes of this paper, we will simply outline a few important technical points. 
Standard normalization techniques may be used if required, though the nilpotent 
structure usually means they are not usually necessary. (The nilpotents used here have 
been multiplied by an extra i for mathematical convenience, and this would be 
removed in normalization.) The momentum or angular momentum operator p, being 
multivariate, can be multiplied with itself to give p , as can the scalar product o.p, and 
so it is possible to use either p or o.p in the quaternion state vector. If fermions can be 
taken as having the four terms in the state vector in a particular order, say, 

(kE + a p + ij m) ; {kE - ii p + ij m); {-kE + ii p + ij m); {-kE - ii p + ij m) , 

then the order of terms in the equivalent antifermion (with the same spin) is also 
determined as: 

(-kE + ii p + ij m) ; (-kE - ii p + ij m); (kE + ii p + ij m); (kE - ii p + ij m) . 

To specify either a fermion or antifermion, it will often be convenient just to specify 
the first of the four terms, with the presence of the other three automatically 
understood. We will, for example, refer on occasions to 

(kE + ii p + ij m) (-kE + ii p + ij m) 

being a scalar on the basis of this understanding. 

Reversing the fermion or antifermion spin state requires reversing the signs of p in 
each term. Boson state vectors are constructed from the scalar product of a fermion 
and an antifermion state vector, arranged as bra and ket (in either order), for example, 

(kE + ii p + ij m) (-kE + ii p + ij m) + (kE - ii p + ij m) (-kE - ii p + ij m) 

+ (-kE + ii p + ij m) (kE + ii p + ij m) + (-kE - ii p + ij m) (kE - ii p + ij m) , 
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and are always pure scalars (in this case, SE before normalization). The expression, 
here, gives the general case; many bosons, of course, will have m = 0. It is also the 

expression for a vector boson. A scalar boson would have (kE + ii p + ij m) (-kE - 
a p + ij m) ... , and cannot have m = 0, as (kE + ii p) (-kE - ii p) . . . then disappears. 

A vacuum operator applied simultaneously to all four solutions is most 
conveniently represented by a diagonal matrix, premultiplied by a row state vector or 
postmultiplied by a column state vector. The vacuum wavefunction operator (when 
applied to a row vector) is always k x matrix form of state vector X exponential term: 



^E + Hp + ijm 

kE- Hp + ijm 
-kE + Hp + ijm 



V 











-i(Et - p.r) 



-kE - Hp + ijm J 



The order is reversed when applied to a column vector. The vacuum operator omits 
the exponential term. Applying the vacuum operator leaves the form of a fermion or 
antifermion quaternion state vector unchanged. It appears to be significant that the 
infinite series 



(feE + Hp + ijm) k (kE + Hp + ijm) k (kE + Hp + ijm) k (kE + Hp + ijm) ... , 

which may be taken as representative as the individual terms in a fermion state vector 
acting on vacuum, is identical to 

(kE + Hp + ijm) (- kE + Hp + ijm) (kE + Hp + ijm) (- kE + Hp + ijm) ... , 

which may be taken as the creation by a fermion of an infinite series of alternate 
boson and fermion structures, presumably in the vacuum in a kind of supersymmetric 
series, with each particle or particle plus vacuum being its own supersymmetric 
partner. This suggests the existence of a kind of supersjmimetry without a new set of 
particles. 



3 A perturbation expansion of the Dirac equation for QED 



We begin with the Dirac equation for a fermion in the presence of the electromagnetic 
potentials (/>, A. ^ ^ ^ 

k^ + ii a.V + ijm ^= -e (i k</>- i a.A) ^ 



and apply a perturbation expansion to y/, so that 
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where y/Q = (kE + Hep + ijm) e ^'^^ is the solution of the unperturbed equation: 



k^^+ii a.V + ijm 



?«^o = 0, 



and represents zeroth-order coupling, or a free electron of momentum p. Here, of 
course, as always, we take {kE + «o.p + ijm), and the equivalent differential operator, 
as a 4-component vector incorporating the four solutions (± kE ± zi'a.p + ijm) in a 
fixed order, say (kE + «a.p + ijm), (- kE + /io.p + ijm), {kE - iio.p + ijm), {- kE - 
iia.p + ijm). This means that {kE + iia.p + ijm) {- kE + iia.p + ijm) becomes a scalar 
quantity, while {kE + iia.p + ijm) {kE + iia.p + ijm) = 0. 
Using the perturbation expansion, we can write 

k^^ + ii o.V + ijm {y/Q + y/i + y/2 + . . .) = -e {ik^ -i o.A) {y/o + y/i + y/2 + ■ ■ ■) , 



J 



leading to the series 



k-^ + ii a.V + ijm y/o= 



^ dt 



k^ + a o.V + ijm \ y/i= -e {i k^- i O.A) y/Q 



^ a 



k^^ + ii o.V + ijm y/i = -e {i k^- i o.A) y/\ 



Expanding (/ k(l)-i O.A) as a Fourier series, and summing over k, we obtain 



(/ k(^ - i O.A) = E (/ k4) (k) - i o.A (k)) 



so that 



k^^ + ii o.V + ijm \y/^i= -e L (/ k^ (k) - i o.A (k)) e' y/o 



= -eZ {i k^{k) - i o.A (k)) e'"^" {kE + iia.p + ijm) g-'^^'-P-'-) 
= -eL {i k^ (k) - i o.A (k)) {kE + ao.p + ijm) e''^^' ' * ^^''^ . 



Suppose we expand y/^ 



y/i=I. vi{E, p + k)e 



-i(Et - (p + k).r) 



Then 



kf^+ii o.V + ijm v,{E, p + k) g-'^^'" * ^ 



= -ei:{ik(l> (k) - i o.A (k)) {kE + Ha. p + ijm) e^^' " ^ '^^-''^ 



L (A:^ + «o.(p + k) + ijm) vi(£', p + k) e 



A{Et - (p + k).r) 



= E (/ k(t>(yC) - i O.A (k)) {kE + «o.p + i/m) ^-'(^'-(p + 
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and, equating individual tenns, 

{kE + Ha(p + k) + ijm) v\{E, p + k) = -e {i A;<A(k) - / o.A (k)) {kE + j/ap + ijm) . 
We can write this in the form 

vi(E, p + k) = -e [kE + ii a.(p + k) + ijmY^ (i k(p (k) - i o.A (k)) (kE + iio.p + (/m) , 
which means that 



This is the wavefunction for first-order coupling, with an electron (for example) 
absorbing or emitting a photon of momentum k. Applying this to the second term in 
the perturbation series, we obtain: 



L (-e) [kE + ii a (p + k) + ijmT\i kcj) (k) - / o.A (k)) {kE + ao.p + ijm) e''^'^' " * ^ "^■''^ . 

We now introduce Z {i A:^(k') - i o.A (k')) e^''^ for the 4-potential involved with 
this term, and summing over k and k'. 



y/i = -el. [kE + ii a.(p + k) + ijm] (i k^ (k) - / o.A (k)) 



X (kE + «a.p + i/m) e 



-i(a-(p + k).r) 




^«/^2 = Z L V2(£:, p + k + k': 



')e 



-i(a-(p + k + k').r) 



and obtain 




-i(£f-(p + k + k').r) 



- e L L (j A;<Z> (k') - i a.A (k')) [kE + a a.(p + k) + ijm] 



1-1 



X (-e) (i k<p(k)- i o.A (k)) (kE + /lo.p + (/m) e 



-i(a-(p + k+k').r) 



Comparing coefficients, we find that 



(kE + ii a.(p + k + k') + ijm) V2(E, p + k + k') = -e(ik^ (k') - / a.A (k')) 



X [kE + ii o.(p + k) + ijm] (-e) (i kcj) (k) - / o.A (k)) (kE + HO.p + ijm) . 
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Hence 

V2{E, p + k + k') = + a a(p + k + k') + ijmY^ {-e) {i A:^(k') - / a.A (k')) 
X {kE + a a(p + k) + ijmY^ {-e) (i k(j) (k) - i a.A (k)) {kE + aa.p + ijm) , 

and 

y/2 = 'LL[kE+ a o.(p + k + k') + ijmT^ (-e) (i k^ (k') - i o.A (k')) 
X [kE + a a.(p + k) + ijmT^ (-e) (i k</> (k) - / a. A (k)) 
X (kE + iic.p + ijm) e-'(^^-(P + '^+k').r) 

This is the wavefunction for second-order coupling, with an electron (for example) 
absorbing and / or emitting two photons of momenta k and k'. 

4 Integral solutions of the Dirac equation 

We may suppose that the equation 



V 



k^^ + a a.V + ijm y/i = -e (i k^- i o.A) y/Q 
has an integral solution of the form 

oo 

y/i{x) = / Gi{x, x) y/oix) dx , 

— oo 

where x and x' are 4- vectors, and 

1 °° 

Gi(x, x) = 7rT4 fctp [kE + a o.p + ijm] ^ (-e) (i k^ (x) - i O.A (x)) exp -ip(x- x) , 

—oo 

where p is also a 4-vector. Defining the Fourier transform of (/ k(/>(x')- i a. A (x)) as 

^ oo 

fctp '(ik^(p')-i O.A (p )) exp -ip'ix-x) , 

—oo 

we obtain ^ 

Gi(x, x') = ^2^8 J" J" d^p ctp' [kE + ji o.p + ijmY^ 

X (-e) (/ k(/>(p')-i O.A (p )) exp -i(p +p')(x-x') . 

The procedure can be extended to y/iix) and 

Giix", x') = '^^ff d^p"(^p"' VkE"+ a ap" + ijmY^ 

X i-e) iik(p(p")-io.A(p")) exp -i(p" + p'") (x'-x') . 



Reversing the order of x and x'in the first term, the product of G2G1 becomes 

G2Gi=^^ j-j-j-j-cfpd^p'ctp"ctp"'[kE"+ a o.p" + ijmT^ 

X i-e) (i k(/> (p ") - i a.A (p '")) {kE + ii o.p + ijmY^ 

X {-e) (ik^ip^-i O.A (p )) exp -i(p " + p '") {x "- x ) exp -i{p +p')(x'- x). 

Here, [kE" + ii o.p" + ijmT^ (-e) (i k^Q?'") - i o.A (p'")) [kE + // o.p + ijmY^ {-e) 
(i k</>(p')-i a.A (pj) has the same form as V2(E, p + k + k'). 

5 Renormalization 

In standard QED, renormalization is not needed at the tree level, where there are no 
loops, but is needed when loops are required in the calculation. There are essentially 
three tj^es of graph which remain divergent after the application of gauge invariance 
and Furry' s theorem. These are the electron self-energy graph (with divergence D = 
1), the electron-photon vertex graph (with D = 0), and the photon vacuum polarization 
graph (with D = 2, which can be reduced to D = 0, by gauge invariance). Let us take 
the electron self-energy graph, for a free electron not interacting with any external 
field. In effect, this is a representation of an electron emitting a virtual photon and 
then reabsorbing it. The perturbation expansion for a first-order coupling of this kind 
(for either emission or absorption) produces a wavefunction of the form 

^1 = -e L [JtE + ii o. (p -H k) ijmr^ (ik^ - i o.A) (kE + iio. p + ijm) e''^^' ' * ^ ''^•'"^ , 

where k is the extra momentum acquired by the electron from the photon. If we 
observe the process in the rest frame of the electron and eliminate any external source 
of potential, then k = 0, and (ik(/> - i a.A) reduces to the static value, ik(/> . In this case, 
^1 becomes 

^1 = -e L [kE + ii o.p -1- ijmT^ ik(/) (kE + i/o.p + ijm) g-'(^'-(P+ •'>•'■) . 
Writing this as 

^1 = -e L [kE + ii o.p + ijm]~^ [-kE + ii o.p + ijmY^ (-kE + ii o.p + ijm) ik^ 

X (kE + HO.p + ijm) e"'^^'- * ^ "^^'O , 

we obtain ^1 = 0, for any fixed value of </). In other words, a non-interacting electron 
requires no renormalization as a result of its self-energy. We can extend this idea to 
the photon vacuum polarization, which creates and then annihilates an electron- 
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positron pair. These processes are equivalent, in effect to producing a diagram in 
which an electron and positron travel in mutual opposite directions from or towards a 
vertex with an incoming or outgoing photon, which, in the absence of any external 
potential, is, in principle, the same as an electron emitting and absorbing a photon and 
being undeviated by any momentum change (k = 0), as in the self-energy diagram. 
The calculation of = for the self-energy should, therefore, also apply in this case. 
The electron-photon vertex graph diverges only beyond the tree level, when a photon 
loop links the electron on each side of a vertex formed with an incoming photon. In 
the context of k = 0, this reduces to another equivalent of electron self-energy or 
vacuum polarization. 

It appears, from this calculation that [kE + ii a.p + ijrn\~^ is equivalent to the 
vacuum 'image' of {kE + // o.p + ijm), that is, {-kE + ii o.p + ijm), and this seems to 
be its physical meaning. In the nilpotent formulation, a fermion with quaternion state 
vector {kE + ii o.p + ijm) generates a vacuum image of the form {-kE + ii o.p + ijm) 
by acting on the vacuum operator k{-kE + ii o.p + ijm). However, {kE + ii o.p + ijm) 
k {kE + ii o.p -I- ijm) is, in principle identical to {kE + ii o.p -i- ijm), as is {kE + ii o.p + 
ijm) k {kE + ii o.p -i- ijm) k {kE + // o.p -i- ijm) k {kE + ii o.p -i- ijm) . . ., in an infinite 
series of actions upon the vacuum; and this is mathematically and physically 
equivalent to {kE + ii o.p + ijm) {-kE + ii o.p + ijm) {kE + // o.p + ijm) {-kE + ii o.p 
+ ijm) which generates an infinite series of alternate fermion and boson states. We 
can, therefore, conceive of this in terms of a fermion generating an alternate series of 
boson and fermion loops, which sum up to return to the pure fermion state itself. It is 
effectively a form of supersymmetry in which the supersymmetric partners are not 
new particles but merely vacuum images or couplings of the original state. 

In addition to this, as we have shown in earlier work, the nilpotent representation 
of the Dirac equation is already second quantized because the operator and state 
vector are necessarily identical, and each is quantized. This means that the quantum 
field integrals reduce in principle to the form of a single quaternion state vector.^ This 
representation says nothing directly about the energy value associated with the 
coupling of a charge to an electromagnetic field, but it does define the fact that such a 
single representation must be possible, and, that the infinite number of possible 
fluctuations must be constrained in such a way that they sum up to a single finite 
value for the coupling in any given state. From this perspective, the existence of a 
method for eliminating apparent infinities is no more surprising than the fact that 
complicated physical systems act in such a way that they conserve energy or that the 
complicated mess of quarks and gluons inside a proton somehow 'conspires' to 
deliver a total spin value of Vi to the composite particle. It is also not surprising that 
there is a maximum value of energy involved which ensures that the summation is 
finite. The existence of the quaternion state vector as a single representation of the 
entire quantum field is a statement that a finite summation exists. The calculation at 
the beginning of this section is an illustration of the fact. There is even a mechanism 
for ensuring that it happens in the interactions of quaternion state vectors with the 
vacuum. 
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When we wish to deteraiine the relative strength of the coupling involved in an 
electron interacting with an external field, for example in the process of electron- 
muon scattering or the coupling of the electron to a vector potential to determine its 
magnetic moment, we then need to sum a perturbative series of real interacting terms. 
The process, however, will not be divergent, and it will naturally cut off at the value 
of energy described as the Planck mass (Mp), the energy at which the effects of what 
is usually referred to as quantum gravity take place. This is, we believe, because 
gravity is the carrier of the wavefunction correlations involved in nonlocality, and the 
Planck mass is the quantum of the inertial interactions which result from the effect of 
gravity on the time-delayed nature of the nongravitational interactions, which, in turn, 
produces the inertial mass associated with 'charged' particles. It is, of course, the fact 
that such particles have inertial mass which creates the time delay in the transmission 
of energy; and the inertial mass may also be seen as the result of a coupling to the 
Higgs field, which produces the filled vacuum which allows the instantaneous 
transmission of gravity and wavefunction correlations. Cutting off at the Planck mass 
((he I G)^'^), which is c-dependent is the same thing, in effect, as saying that the real 
interactions (and inertial reactions) are time-delayed or that the fermions involved 
have inertial mass, and it is significant that, in QED, the process actually determines 
the mass-value of the electron. The process is, in effect, creating a c-related 'event 
horizon' such as we would expect for an inertial process as described. 

The significant result from current QED calculations is that the first-order 
correction term for the electric fine structure constant, for purely leptonic interactions, 
at energy of interaction ^, is (1 / 3;r) In (Mp^ I jl), which has the same form as the 
term which we have derived for the quark theory, using the idea of lepton-type 
quarks, and also the correction term which is needed to achieve Grand Unification at 
the Planck Mass. The reasons for the Planck mass being the energy for a l](S) 
unification, involving gravity, have been discussed in earlier papers, together with 
relevant predictions.^'^ 

6 Green's function solution 

To apply the method of Green's functions, we solve for a unit source and then sum 
over the whole distribution. Beginning with 



or 



k^i^ + ii O.V + ijm y/- -e (i k^-i o.A) y/ , 




,ip.x 



where p and x are 4-vectors, we solve for a unit source using 




(x-x") , 
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where Gp is a wave produced at x by a unit source at x'. Translational invariance 
shows that Gp is a function of (x - x'). Then the solution for the whole distribution 
becomes 

y/(x) = -e f ctx' Gfix - x) {i k(/>- i o.A) yKx) 

or 



(kE + «o.p + ijm)e'^'' - e f ctx' Gf(x -x') (ik(/)- i o.A) {kE + i/o.p + ijm)e 



ip.x 



Taking the Fourier transform into momentum space, we obtain 

>.(x-x') 



Gp{x-x') = -^fcfp SAP) e-'P-' 



Then, using the Fourier representation of the ^function, we find 

-ip.{x-x') _ — 1_ f A 



k'^ + ii a.V + ijm Spip) e " " ' = ^^^4 jape 



ip.(x-x') 



from which 
and 



(fcE + /ia.p + i/m) S'iKp) = 1 



1 



^"^^-(kE+iiG.p + ij 



ijm) 



which is the electron propagator in the Fejmman formalism. Conventional theory 
assumes that 



1 ff + m 
p-m p -m 



2 2 

and that there is a singularity or 'pole' (po) where p - m = 0, the 'pole' being the 
origin of positron states. Effectively, on either side of the pole we have positive 
energy states moving forwards in time, and negative energy states moving backwards 
in time, the terms + m) and (- p + m) being used to project out, respectively, the 

2 2 

positive and negative energy states. We then add the infinitesimal term ieXo p - m , 
so that iSpip) becomes 

i{p + m) 
'^'^^^ p'-m' + ie ' 



(jf + m) '^ 
2po 



1 



1 



[Po - ^Ip^ + 171" + ie Po + "Sjp" + 171" - ie 
and take a contour integral over the complex variable to give the solution 

Sp (x -x') = f cfp ^ [-ie(t - n I ^'(x) ^(x') 



+ 



r^. — 2 



+ ie(t'-t) Z ^(x)^(xO] 

r=3 
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with summations over the up and down spin states. 

In the nilpotent formalism, however, the i£ term is unnecessary and there is no 
infrared divergence at the pole, because the denominator of the propagator term is a 
positive nonzero scalar. We write 

1 

Slip) - ^ ^Yo.p + ijm) ' 

and are free to choose 

1 {-kE + u'o.p + ijm) 

(kE + aa.p + ijm) ~ {kE + ao.p + ijm) {-kE + «a.p + ijm) 

_ {-kE + iio.p + ijm) 
" +p^ + mS ' 

which is finite at all values. We can also, for example, write 



J -kE + /fg.p + ijm 



r 1 1 ^ 



{kE + j/a.p + ijm) -2kE [{kE + iic.p + ijm) { -kE + iic.p + ijm) 

and other, similar, forms. Our integral is now simply 

SF{x-x') = jd'p-^^0{t-t')^{x)^{x'), 

in which 

^{x) - {{kE + a p + ij m) ...) exp {ipx) , 

where {{kE + it p + ij m) . . .) is the bra matrix with the terms: 

{kE + a p + ij m); {kE - ii p + ij m); {-kE + ii p + ij m); {-kE - ii p + ij m) , 

and the adjoint term becomes 

^(x') = {{kE - ii p- ij m) ...) {ik) exp {-ipx) , 

with {{kE - ii p - ij m) ...) {ik) a ket. The reason for this success is apparent. The 
nilpotent formulation is automatically second quantized and the negative energy states 
appear as components of the nilpotent wavefunction on the same basis as the positive 
energy states. No averaging over spin states or 'interpreting' -E as a reversed time 
state is necessary; the 'reversed time' state occurs with the t in the operator ^, and 
there is no need to separate out the states on opposite sides of the pole. 

We can use the electron propagator to define the photon propagator. 
Conventionally, we derive the photon propagator, Af{x - x'), directly from the Klein- 
Gordon equation, while recognizing that its mathematical form depends on the choice 
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of gauge. In the presence of a source field, represented by current j(x), we can write 
the Klein-Gordon equation in the form 



" r72 2 



^x) =j(x) 



and, with Af{x - x') as Green's function, we have 



Af(x-x')^-S'^^\x-x) , 



for which the solution is 



(/)(x) = (/Jb(x) - J (fx' Af (x - x') j(x') , 



where (/)o(x) is a solution of the equation in the absence of sources: 



dt 



-^,2- y -m 



^x) = . 



Using the Fourier transform 



Af (x - xO = f ctp exp -ip{x - x) Af (p) 



and applying the operator 



df 



-m 



to both sides, we obtain: 



J 

Af{x - X') = 2 2 

^ p -m 



Clearly, the relationship of the electron and photon propagators is of the form 



or, m our notation. 



Sf (x - xO = (/ fd/u + m) Af (x - xO , 



Sf (x - x') = ((kE + a p + ij m) .. .) Af (x - x') 



This is exactly what we would expect in transferring from boson (Klein-Gordon field) 
to fermion (Dirac field), using our single vector operator. Now, using 



iSpip) = 



1 



1 



^Po [pq - ^Jp^ + rn' + ie po + ^jp^ + - ie 

which is the same as the conventional electron propagator up to a factor + m), we 
can perform virtually the same contour integral as in the case of the electron to 
produce ^ ^ 



1 



IAf (x - x') = / d'p ^^(t- n <Kx)f(x') 
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where CO takes the place of E I m. This time, of course, <f(x) and are scalar 
wavefunctions. In our notation, they are each scalar products of the 4-component bra 
term {(kE + ii p + ij m) ...) and the 4-component ket term ((-kE + ii p + ij m) ...), 
multiplied respectively by exponentials exp (ipx) and exp (ipx), expressed in terms of 
the 4-vectors p, x and x'. In the nilpotent formulation, </)(x)</f^{x') reduces to a product 
of a scalar term (which can be removed by normalization) and exp ip(x -x'). 

Now, the product ^(x) ^(x'), for the electron propagator, can be interpreted as the 
product of ((kE + ii p + j/ m) . . .) exp (ipx) and the 'vacuum operator', ik ((kE + ii p + 
ij m) ...) exp (-ipx'), which reduces to the product of a scalar (again removed by 
normalization), ((kE + ii p + ij m) ...) and exp ip(x - x). The integrals, like the 
propagators, are related by the factor ((kE + ii p + ij m) ...), which suggests that we 
should define the photon propagator, in terms of the electron propagator, for which, in 
the nilpotent formulation, we can do a completely non-divergent integral without 
additional infinitesimal terms, rather than as a result of the Klein-Gordon equation. In 
principle, this must reflect the fact that photon processes cannot ultimately be 
described independently of processes involving paired fermions or fermions and 
antifermions (one of which may be manifested indirectly in terms of a potential, rather 
than directly as a charge). We can then write 

(kE + ii p + ij m) Ap (x-x') = Sp (x-x') = (-kE + ii p + ij m) , 

—2 

which means that (after normalization) Sp (x - x') takes the expected form of ^ = 1 / 

2 2 2 

(E - p ), with m = 0. The apparently divergent nature of the photon propagator 
integral, in conventional theory, is simply a reflection of the fact that photons, unlike 
electrons, have no independent existence and are not conserved objects. However, if 
we define the photon propagator in terms of the nilpotent fermion propagator integral, 
we will be using the automatic cancellation process which the nilpotent representation 
of the fermion introduces. 

7 The propagator method in lowest order 

The Dirac equation in the presence of external Coulomb potentials, 

( 5 ^ 

k^+ii a.V + ijm \f/= -e (i k<^-i o.A) y/ 
may be solved by introducing the Dirac propagator, so that 

^ a ^ 

V 



k^^ + ii a.V + ijm Sp (x-x') = -e (i k^- i a.A) S^'*^^ (x - x) 



Then 



y/(x) - \i/q(x) + J (fx' Sp (x - x') (-e) (i k(/>- i o.A) y/. 
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where y/dix) is the solution of 



k^^+ii a.V + ijm U^b = 



The Dirac equation with potentials may also be written in terms of the Hamiltonian, 



H = ik ii a.(V + i eA) + ijm = Ho + Hi , 



where 



and Hi is the interaction term. 



Then 



^0 = i i tJ'V + ijm 



Hi = ik e(/>- i e O.A) 



k^^-^iH 



and 



3 ^ 
k^^ + iHo + i Hij G(x -x') = S^^'^ (x - x) 



Solving the Green's function, G(x - x"), for the general case involving interactions, 
allows us to use Huygens' principle to find the time evolution of the wave, by 
summing all the secondary wavelets produced by infinitesimal point sources on the 
wavefront. In mathematical terms (changing over the variables from 4-vector 
notation), 

V<x, = J (f'x'Gfix -x',t- O yKx\ O . 



Using the symbolic identity 



G = Go + GqHiGo + GqHiGoHiGo + 



and a power expansion, and applying Huygens' principle, the time evolution of the 
wavefunction becomes described by 

^Kx, = y/o(x, t) + j- AiGo(x -xi,t- ti) Hi(xi, h) y/oO^i, h) 

+ J cfxi cfx2 Go(X -Xi,t- h) Hi(Xi, ti) Go(Xi - X2, h - t2) y/Q{X2, ?2) + ... . 

To calculate the transition probability, in lowest order, that a free plane wave in initial 
state / emerges in final state /after scattering off a potential, we multiply the equation 

yAx, t) = (jhix, r) + J A' Go(x -x\t-t') Hi{x\ t) (/h(x', t) 
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on the left by ^* and integrate, obtaining the scattering matrix, 

Sfi=Sri + i J(fx ^(x) Hiix) ^(x) + . . . 

Considering, now, an electron interacting with external Coulomb potentials, we 
introduce 

y/(x) = y/o(x) + J cfx' Sf (x - x') (-e) (i kcj)- i o.A) y/, 

and obtain 

yX^x) = y/lx) + (-e) J cfx' 0(t- t')j (fp -^^^ ^ yKx) yKx')(i k^x) - i cAixy) yAx') . 

Multiplying both sides of the equation from the left by y/f and integrating over all 
space and time, we obtain the scattering matrix, to lowest order: 

Sfi = Sfi -ie S y/f y/i ctx . 

8 Electron scattering 

From the previous section, we find that the amplitude for the scattering of an electron 
from initial state y/i to final state ^is given by 

Tfi = ieS y/f YfAn Wi d^x 

- ie f ik (kEf+ j/o.p/+ ijiUf) exp (/ (Ejt - p/.r)) 

X (ik^ - i a. A) (kEi + iic.pi + ijirii) exp (/ (Eit - Pi.r)) d rdt 

= ie f ik (kEf+ iicpf+ ijrrif) (ik^ - i a. A) 

X (kEi + iic.pi + ijmi) exp (- i ((£■,• - Ef)t - (p, - p/).r)) d rdt . 

It is convenient to write this in terms of the electromagnetic transition current 
between the initial and final states, 

but this is only possible if we redefine the scalar product J0'"as a product of diagonal 
matrix (y^) and column vector (A""), to separate out the terms ik and -io, within j^, 
which have quite different actions on (kEi + «o.p/ + z/m,). In effect, k reverses the sign 
of E (time-reversal) and i reverses the sign of p (parity), so transforming into a 
matrix is like applying the metric tensor, g^^. Thus, the row vector now becomes a 
diagonal matrix of the form: 



16 











A 





-i(7i 














-io2 













-iC73 J 



It is convenient, also, to write exp (- i ((Ei - Ef)t - (p, - P/).r)) = exp (-i qjc), in 
tenns of a 4-vector difference x and a 4-momentum difference q, of which the vector 
part is q = p, - p/. So, the electron current becomes 

j^ = -e y/f y/i = - e ik (kEf+ iicpf+ ijrrif) yju (kEi + ua.p,- + //m,) e~ "''^ . 

Defining A\x) as the 4-vector potential associated with the static charge, and A^(q) as 
its Fourier transform, with vector component A'"(q), 

A\q)= ^ xe- A\x)(tx , 

the amplitude for the process becomes 

Tfi = ie fy/f y/i A^(x) A = - / fif^^^x) A , 

= - ek (kEf+ iia.pf+ ijmf) {kEi + HO.p,- + i/m,) A%q) . 

For a static source, with time-independent A'"(;c), 

A%q) = J exp (- / {Ei - Ef) t) dt J e'*'' A^(x) A . 

= l7td{Ef-Ei)A''{fi) . 

Applying Maxwell's equations for time-independent A'"(x), we have 

VV(x)=-/(x) , 

and so 

J e'l-'' A\x) A = J (V' A^(x)) e'l-" A 
= -|qpA^(q) = -/(q) . 

Hence 

A^(q)=[^/(q) > 

and 

r^- = - 2i7rS{Ef- Ei) e k {kEf+ iia.pf+ ijmf) y^ {kEi + iia.pi + ijmd ri]2/(q) • 
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Removing the ^function, we obtain the co variant amplitude: 

-i M = -ek (kEf+ iiG.pf+ ijmf) y^i (kEi + iia.pi + ijnii) j"^7'"(q) . 

Applying energy conservation, Ef= Et, and q' = -\ qP, which means that the invariant 
amplitude may be written: 

-iM = -ek (kEf+ «a.p/+ ijnif) f{kEi + jio.p, + ijmi) (-^/^(q)) • 

We can, of course, use this equation, for example, with = )f) and g^y = 1, and iZe 
replacing ij (q), to find a value for the angular distribution (|M1 ) in the case of 
Rutherford scattering off a static nuclear charge Ze; and we can also assume that more 
complicated calculations can be done using the Feynman rules which codify the 
perturbative method. However, our main concern is with developing an approach to 
renormalization. Let us suppose, therefore, that the exchanged photon fluctuates into 
an electron-positron pair. Application of the Feynman rules requires a factor (-1)" for 
a diagram containing n fermion loops (in this case, n = 1). It also requires an integral 
over ctp I {iTtf to sum over all possible values of the unobservable p, in the part of 
the invariant amplitude referring to the loop. The invariant amplitude is, therefore, 
now: 

-/ M = -(-)^ e k (kEf+ iicpf+ ijmf) fikEi + HO.p, + i/m,) 

1 n\ i{kE+ii<3.\t + iim)Bx iiiio.jip - (\) + ijm) bx a 

^ / ^''^^"P (E' - p^-J) ((p - qf -mS 

X=¥^(-r(q)). 

In principle, the addition of this term for the photon loop to the lowest-order invariant 
amplitude for electron scattering can be achieved by the addition of a modifying term 
to the lowest-order propagator, so that 

q q q q 



with 



q q ^ q 



V (9 ) = (-) ^ / Tr \{ier,) ^E'-p'-m') ^''^^^ {ip-q)'-m') ) 

i6 



Bjorken and Drell and others show that, with the omission of terms which 
disappear when the propagator is coupled to external charges or currents, I^y can be 
written as 

Inv = -gnvq^Kq^) 
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with 



CXD 1 



1- m' 1^^' 



m 

where m is the electron mass. In principle, the first term of this integral is divergent, 
but with Mp representing a natural cut-off value for the mass, I(q^) becomes the 
convergent: 



Mp 



1 



2, ^ f dp 2a I 



1 - 1 
m 



m 







For small values of (-q ), we have 



Kq)-^ In — r + TT- ^ 
^ in tn IStt m 



while large values of {-q ) lead to 

In 



^ 2^ « , Mp^ a-q' a , Mp^ 
% ) 



T T-'T- ~% = ~ In T 

2>n m Stt m Stt -q 



The first term in each case may be considered as a modification to the fine structure 
constant or the electric charge value as produced by the lowest-order Feynman 
diagram. The modified (or measured) charge (e') then becomes: 



e - e 



e ivip 

l-T' In — 1 
y 37t m J 



1/2 



In principle, of course, 'charge' really represents an existence condition rather than a 
value. The 'value' of charge is a measure of the coupling to the energy of the field. 
Charge, as a fundamentally imaginary quantity, has no value. Its 'value' can only be 
determined by its coupling to another charge, and then is only determined in terms of 
the energy of the coupling. In this respect, charge is unlike mass, which, as a real 
quantity, has a value independent of its coupling with other masses. We can only 
know that a charge exists by its interaction with other charges, but we can know of the 
existence of a mass through its inertia, irrespective of the existence of other masses. 
The scaling of the interaction value associated with an electromagnetic charge at 
different energies is, therefore, an expected process, which is related to a similar 
scaling of the electron mass (calculated from the self-energy diagram), and stems 
from the same vacuum origin by which this is created. In the absence of an ultraviolet 
divergence in the integral, it produces the expected result, demanded by the nilpotent 
algebra, of a finite (but not independently determined) value of charge at all energies. 

It is significant that the alternative (BPHZ) method of renormalization, via the 
adding of counterterms to the action, results in the same scaling of mass and charge as 
the direct multiplicative method described above. We can consider this as equivalent 
to the automatic addition of inherent supers ymmetric partners in the nilpotent algebra. 
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which has the advantage of being a required property at each successive order. A 
similar process of cancelling the contributions from electrons and those from 
electrons plus very soft (undetectable) photons has been used to remove the infrared 
divergences in bremsstrahlung and pair production. Such procedures can be related to 
our analyses of divergent diagrams in section 5, and our discussion of the creation of 
alternate fermionic and bosonic vacuum 'images' in section 2 (where the additional 
photon energy is actually zero, and so necessarily undetectable). 



9 Strong and weak analogues 



Both strong and weak interactions have relationships between fermion and boson 
propagators which follow the same overall pattern as that in QED. The quark 
propagator in the strong interaction is simply that for any fermion, and the fermion 
propagator, Sf(p) or iSpip), applies to the weak, as well as to the electromagnetic 
interaction. The conventional 

•c f ^ ' '(^ + 

p-m p -m 



becomes 



i(-kE + g'o.p + ijm) 
(E' +p' + rn) 



in the nilpotent formulation, and its general applicability means that we can eliminate 
the infrared divergence in the fermion propagators for all of the interactions. 

It is interesting that the 'renormalizability' of the combined electroweak interaction 
is related to the very mechanism which gives masses to the fermions and gauge 
bosons. For obvious reasons, a quantum field integral taken over all values of p will 
only be finite, as the nilpotent algebra demands, if the index of p in the integrand (or 
divergence D) is less than 0. Now, the propagator for the combined electroweak 
gauge bosons is 



1 



m 



-g/zv+(l-^ 



p - 



JJ 



where ^ is the 't Hooft gauge term, which appears in the gauge fixing term in the 
Lagrangian for the interaction: 

1 



This term removes the unphysical (massless scalar) Goldstone boson ^, which arises 
from the spontaneous symmetry breaking produced by the filled weak vacuum used to 
eliminate negative energy states. If ^ is finite, then as oo, A^^ p~^, like the 

pure photon propagator, which, in the absence of any gauge choice, becomes: 



JJ 
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However, for ^ — > oo, we have the propagator for a massive vector boson theory 
without massless component. 



1 



p -m 



-guv 



PuPv 
2 

m 



which becomes a constant when — > oo, leading to infinite sums in the diagrams 
equivalent to those in QED. One of the most convenient choices of gauge is ^ = 1 
(Feynman gauge), which leads to an electroweak boson propagator. 



p -m 



entirely analagous to that for the photon in the same gauge, and similarly linked by 
the factor {kE + zip + ijm) to the fermion propagator, as in QED. (It may be possible, 
here, to link the existence of massive weak bosons to the creation, in the nilpotent 
representation, of massive bosonic states via the interactions of fermions with the 
vacuum.) 

The gluon propagator in QCD is again of the same form as that of the other gauge 
bosons: 



-guv 



which reduces, in Feynman gauge, to 



and which only differs from the electroweak boson propagator by a factor, 3^^, arising 
from the fact that freely propagating gluons have fixed colours. 



10 Conclusion 



It is not necessary to follow in any further detail the procedures of QED, or derive 
further results by examining specific Feynman diagrams. The outline derivations in 
this paper show that the nilpotent method is not only amenable to QED calculations, 
but also has something additional to offer. It removes the infrared divergence from the 
integral for the fermion propagator by, in effect, combining the fermion and 
antfermion solutions in a single expression. By being already second quantized it 
requires the existence of finite values for fundamental quantities such as mass and 
charge at all energies, whatever integrals are performed to incorporate vacuum 
fluctuations, and the algebra itself provides successful methods for the automatic 
removal of all possible divergences. Since, the nilpotent formulation does not 
incorporate direct information about the measure of the coupling of charge to the field 
energy, then there is no fixed 'value' of charge to be renormalized, and perturbation 
expansions merely correlate the different measures of coupling which occur at 
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different external field energies. Even when perturbation expansions need to be used, 
the theory removes divergences by demanding that there must be a finite upper limit 
to the energy used, which other work on the algebra tells us must be the Planck 
mass. ' Because the process involved (via a filled vacuum) is the one which actually 
creates the masses of fermions, then our measure of the interaction of a charge with 
the field energy becomes simply an expected scaling mechanism, as happens also in 
classical contexts, rather than a desperate means of avoiding otherwise catastrophic 
results. It would appear that both the concept of 'renormalization' and the problem of 
divergences may be eliminated by using the nilpotent algebra, though the successful 
mathematical structure of QED remains intact. 
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